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ABSTRACT 

Context. Microwave (MW) and hard X-ray (HXR) data are thought to be powerful means for investigating the 
mechanisms of particle acceleration and precipitation in solar flares reflecting different aspects of electrons interaction 
the ambient particles in a presence of magnetic field. Simultaneous simulation of HXR and MW emission with the same 
populations of electrons is still a great challenge for interpretation of observations in real events. Recent progress in 
simulations of particle kinetics with time-dependent Fokker-Planck (FP) approach offers an opportunity to produce 
such the interpretation. 

Aims. In this paper we apply the FP kinetic model of precipitation of electron beam with energy range from 12 keV 
to 1.2 MeV to the interpretation of X-ray and microwave emissions observed in the flare of 10 March 2001. 
Methods. The theoretical HXR and MW emissions were calculated by using the distribution functions of electron 
beams found by solving time-dependent Fokker-Planck approach in a converging magnetic field (Zharkova at al., 2010; 
Kuznetsov and Zharkova, 2010) for anisotropic scattering of beam electrons on the ambient particles in Coloumb 
collisions and Ohmic losses. 

Results. The simultaneous observed HXR photon spectra and frequency distribution of MW emission and polarization 
were fit by those simulated from FP models which include the effects of electric field induced by beam electrons and 
precipitation into a converging magnetic loop. Magnetic field strengths in the footpoints on the photosphere were 
updated with newly calibrated SOHO/MDI data. The observed HXR energy spectrum above 10 keV is shown to be 
a double power law which was fit precisely by the photon HXR spectrum simulated for the model including the self- 
induced electric field but without magnetic convergence. The MW emission simulated for different models of electron 
precipitation revealed a better fit (above 90% confidence level) to the observed distribution at higher frequencies for 
the models combining collisions and electric field effects with a moderate magnetic field convergence of 2. The MW 
simulations were able to reproduce closely the main features of the MW emission observed at higher frequencies: 
the spectral index, the frequency of peak intensity and the frequency of the MW polarisation reversal while at lower 
frequencies the simulated MW intensities are lower than the observed ones. 

Key words. Sun: fiares — Sun: X-rays, gamma-rays — Sun: radio radiation — scattering — polarization 



1. Introduction jKundu' '1985': 'Kundu et al." '2001a"b', '2004'; 'Vi lmer et all 

^ [2002; Wilson & Holnian 2003; Kundu ct al. 200^. 

^ Spatial configurations of flaring sites in hard X-ray (HXR) 

- ■ ■' and microwave (MW) emissions are highly variable, they However, there are substantial discrepancies in the loca- 

may be formed by a single loop with one coronal source and tions of these HXR and MW sources observed in the same 

two footpoints ((Masuda et al.l 11994 iKundu et all l2001b[ ) Aare — e.g., they are often separated from each other, and 

or sever al sets of loops with a few footpoin ts and coronal the areas covered by each emission are substantially dif- 

sources (|Hanaokal 119961 : iKundu et alll2001al) leading to dif- ferent. Observations often suggest that the sources of MW 

ferent mo dels for d ifferent types of flaring events (see the and HXR emission may be separ ated in depth s within a 

review bv iMelros^ (jiggQ) ). flaring loop (Kundu 1985; Takak ura et al.|[l995l : ISui et all 

I2002D . The intensities of HXR and MW emissions in the 

Simultaneous observations of hard X-ray emission opposite legs of the same flaring loop can be strongly anti- 

and microwaves in footpoints of solar flares often show correlated being in one footpoint of the same loop higher 

their close temporal correlation pointin g out to thei r in MW and lower in HXR emission, while in the other 

common origin ( Aschwandcn 2005; Bas tian et al.l Il998l) . foot point be ing highe r in HXR and lower in MW emis- 

There is a high likelihood that these emissions are pro- sion (lKunduili l985: Ta kakura et al.lll995l:lKundu et al.ll2004l: 

duced by the same population of non-thermal electrons iGrechnev et al..>2008i) . In addition, there is sometimes a de- 
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lay b etween HXR and MW em i ssions occurring in the same 
flare (|Kundu et al.lll995ll2004l: ISui et al.ll2002h . 

MW emission demonstrates a high variability of po- 
larization from a few pe r cent (lAltvntsev et alJ l2000t 



20031) UP to 



Lee fc Gary! 



[Fleishman fc Melnikovl [20031: iFleishman et alT 
50-100% iDulk igsa lAltvntsev et alJ 12005 
I2OOOD . Unlike HXR emission, which properties are con- 
trolled by the parameters of emitting electrons and am- 
bient plasma particles and Comp ton backscattering, i.e. 
albedo effects (|Kontar et all I2006D , MW emission is de- 
termined in addition by various radiative transfer ef- 
fects. Nevertheless, general perception of MW polariza- 
tion is that highly collimated electron beams tend to 
prod uce MW emission with a higher polarization de- 
gree fPleishman & MelnikoV '2003; 'Fleishman et al.| l2003t 
[Mclnikov ct al. 2008; Fleishman k Kuznctsov 2010). 

Moreover, the resulting photon spectra of HXR and 
MW emission produced during flares are also significantly 
different. In particular, HXR emission from strongest flares 
often has the elbow-type spectra with double power-law 
energy distributions (Lin et al. 2003; Holman et al. 2003) 
and noticeable differences between the spectral indices at 
lower (below 70 keV) and higher energies. MW distribution 
in frequency ()Dulkl Il985t iBastian et al.l I1998D has a maxi- 
mum with a gradual decrease of the intensity towards lower 
frequencies (optically thick emission) and a negative power 
law distribution towards higher frequencies (optically thin 
emission). The peak frequency usually varies from 3 to 20 
GHz for different flares or different sta ges of the same flare 
(|Nita et al.l[200l iMelnikov et al.l l2008h. 

These differences clearly indicate different transport sce- 
narios for high-energy particles precipitating into the foot- 
points of a flaring loop. The mechanisms of transport affect- 
ing HXR and MW emissions are well known to be substan- 
tially different: MW radiation is related to gyrosynchrotron 
emission of high-energy elect rons with energies f rom few 
tens keV (iKundu et alTl2001a^ up to several MeV (jBastianl 
ll999HKundu et al.ll2004D . while HXR radiation is often pro- 
duced by the electrons with much lower energies from 10 
to 30 keV (see, for example, iLin et alll2003l : iHolman et al.l 
l2003f) . In spite of a large number of simulations for the in- 
terpretation of HXR and MW emi ssion with thick target 
models (see , for example, reviews Fleishman & MelnikoV 
2003; Fleishman et al.ll2003t [Fleish man fc _ Kuznctsov 2Q10e 
Krucker et al. 20081 l2010l and refer ences therein), so far 



there are only a few simul ations (jMelnikov et al.l l2008t 
iKuznetsov fc Zharkoval l2010h considering the problem of 
particle transport in other than collisional energy loss mech- 
anisms and their effects on HXR and MW emission. 

However, recent kinetic models calculated electron 
distribution functions numerically by solving the time- 
dependent Fokker-Planck kinetic equation for precip- 
itation of electron beams into a flaring atmosphere 
with converging magnetic field. Beam electrons with 
an energy range from 12 keV to 10 MeV are as- 
sumed to lose their energy and change propagation di- 
rections in collisions and due to O hmic losses in the 
electric field induced by the beam (iDiakonov fc SomoyI 
1988t IZharkova et al.1 119951: IZharkova fc Gordovskvvl 120061: 



Battaglia fc Ben J 2008 ) and due to scattering in converg- 
ing magneti c field jLeach fc PctrosianI IT98lt iMcClementsI 
[l992a,b; Si verskv fc Zh arkova 2009). These solutions were 
applied for the explanation of HXR emission produced by 
solar fiares which provided a very good fit for HXR in- 



tensities (lMcClementslll992bt ). directivity and polarization 
([Zharkova et al.ll2010t) . Recent extension of these solutions 
to the interpretations of MW emission from solar flares ap- 
pears to pr oduce promising results in c loser fitting to ob- 
servations ()Kuznetsov fc Zharkoval[2010() . 

These kinetic models can be further tested by the si- 
multaneous interpretation of HXR and MW emission ob- 
served in the same fiare. For this purpose we selected 



studied (iLiu et al. 


2001": 'Bind 20031: Uddin et al.l 2004: 


Chandra et al. 12006 


: Altvntscv et al.ll2008l 


Mclnikov et al.] 



2008) because of the following reasons: 1) it had produced 
both kinds, HXR and MW emissions, MW with a pro- 
nounced polarization; 2) it was a rather powerful flare which 
had to be produced by powerful electron beams allowing 
one to explore the effects of various energy loss mechanisms 
on HXR and MW emissions; and 3) it was located far from 
the disk centre that allows us to compare the energy loss 
effects in a flaring atmosphere seen by an observer from the 
Earth at larger viewing angles. 

The first detailed analysis of a spatial dynamics of 
this flare in HXR and MW emissions was made by 
IChandra et al.l ([20061 ). They showed that the spatial struc- 
ture derived from the microwave emission allows one to 
assume an interaction between two loops, one of which was 
a small, newly emerging loop, while the other was a large 
overlying loop, similar to the events identified by Hanaoka 
(iHanaokal [19961 [l999allbl) 

Later lAltvntsev et al.l ([2008fl also interpreted the high- 
frequency MW emission in this event produced by elec- 
trons precipitating into a homogeneous source or into two 
homogeneous sources located in two loop legs where the 
depth variations of the parameters of high energy electrons, 
plasma density and magnetic field were neglected. The au- 
thors though concluded that such electrons were rather 
anisotropic, or beamed, because (i) the microwave emission 
consisted of many short (the order of seconds) broadband 
pulses supposedly marking quite a few elementary acts of 
beam injection into the loop legs, (ii) the microwave emis- 
sion was 0-mode polarized at 17 GHz (optically thin part 
of the spectrum) meaning that the observed MW emission 
was an immediate outcome of the beam precipitation rather 
than the radiative transfer effects, (iii) type Ill-like drift- 
ing bursts were observed at lower frequencies meaning that 
beam electrons streamed to higher atmospheric levels. 

The goal of the current research is to extend the in- 
vestigation of this flare to simultaneous interpretation of 
HXR and MW emissions and to fit their observed energy 
distributions with the simulations of electron precipita- 
tion into a converging magnetic loop by using the Fokker- 
Planck kinetic approach taking into account collisional and 
Ohmic losses in the electric field induced by well collimated 
beam electrons. Observations are described in sectionj^l the 
model and method of simulations are described in sectionjSl 
the results of the model fit to observations are shown in sec- 
tion 21 and the conclusions are drawn in section [51 



2. Observations 

2.1. Instrumentation 

We used microwave total fiux reco rds from N ob eyama 



Radio Polarimeters (NoRP; iTorii et all 119791: 



IShibasaki et al.l 119791: iNakaiima et al.l I1985D at 
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Fig. 1. The 10 March 2001 event. Time profiles of the average brightness temperature at 1-80 GHz (NoRP) in total 
intensity (a) and polarization (b). (V/I)max is the maximum degree of polarization at the central peak (04:03:40). The 
numbers in Y axes are in the solar flux units (sfu). 



3.75, 9.4, 17, 35, and 80 GHz. The time resolution of 
routine NoRP data available at the NoRP Web site is 
1 s, and for flares it is 0.1 s. Nobeyama Radioheliograph 
(jNakaiima et al.l[l994l) produces images at 17 and 34 GHz. 
The imaging interval of the flare mode data which we used 
in this work was 0.1 s. 

The line-of-sight magnetograms have been produced 
with the Michelson Doppler Imager (MDI; IScherrer et al.l 
Il995f ) on SOHO. The full disk magnetograms with a res- 
olution of 1.98" re-calibrated in December, 2008 were ob- 
tained from the Solar Data Information Centre at Stanford 
University, US. 

Information about HXR emi ssion was supplied by 
instruments on Yohkoh satellite (jKosuei et al.lll991al lbh. 
Images of HXR sources and fluxes with a high temporal res- 
olution (0.5 s) were obtained from Hard X-Ray Telescope 
(HXT) data. HXT carried out observations in four spec- 
tral bands: L, 14-23 keV; Ml, 23-33 keV; M2, 33-53 keV; 
H, 53-93 keV. The spectra were obtained from data of the 
Hard X-ray Spectrometer (HXS), w hich was one sensor of 
the Wide Band S pectrometer (WBS: rYbshimori et al.|[T99lt 
ISato et al]l2006f) . HXS data provided hard X-ray spectra in 



a wide energy range of 20-657 keV. The temporal resolution 
of HXS data was 4 s, i.e., slightly lower than the temporal 
resolution of HXT. 

2.2. HXR and MW bursts and their sources 

The 1B/M6.7 flare of 10 March 2001 associated with a CME 
occurred in active region 9368 (N27W42) and was previ- 
ously studied in detail in a number of papers listed above. 
However, recent re-calibration of magnetic field strengths 
in SOHO/MDI magnetograms alerted a review of previous 
results and inspired an analysis of additional observations, 
which were not yet considered. 

2.2.1. Light curves of MW and HXR emissions 

As found bv lAltvntsev et al1()2008l ). broadband short pulses 
were present in the whole frequency range 1-80 GHz of 
NoRP (see time profiles in Fig.[T^). The MW time profiles 
in total intensity at different frequencies were very similar. 
Three main peaks are detectable in the MW burst with a 
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Fig. 2. Light curves of the HXR burst registered in the 
HXT/M2 channel (33-53 keV) and the HXS record inte- 
grated over the whole band (20-657 keV, histogram) along 
with the microwave burst registered by NoRP at 17 GHz. 
The vertical dotted line marks the main peak under con- 
sideration at 04:03:40 UT. 



total duration of about 40 s (see Fig. [TJi). The first MW 
peak at about 04:03:30.90 {all times hereafter are UT) was 
the lowest one, the second (central) MW peak at 04:03:40 
with a total duration of about 5 s was the highest one. The 
third MW peak had a double structure with two sub-peaks 
occurring at 04:03:44.90 UT and 04:03:46.80 UT, respec- 
tively. 

The time profiles of the polarized MW emission (Stokes 
V component) around the central peak are shown in 
Fig. [Hj. They coincid e with the light curves shown by 
lAltvntsev et all (|2008[ ). The MW emissions of all the three 



peaks had a positive (right-handed) circular polarization 
(RCP) at low frequencies (1-3.75 GHz), and their polar- 
ization at higher frequencies (9.4-35 GHz) was negative 
(left-handed, LCP), as Fig. [T)d shows. 

In order to understand the dynamics of high-energy par- 
ticles in this flare, let us also use the whole energy range of 
the HXR emission recorded by Yohkoh/HXS and estimate 
parameters of the electron beam, th us ex tending studies 
carrie d out by lAltvntsev et al.l (|2008[ ) and IChandra et al.l 
()2006( ) that were confined to the MW emission. 

This fiare was observed by Yohkoh/HXT in four energy 
channels L, Ml, M2, and H. Thus, the available energy 
bands of HXR channels allowed us to derive photon spec- 
tra up to 600 keV. The HXS light curve presenting the 
integrated flux within 20-657 keV range is plotted along 
with HXT and NoRP data as shown in Fig. H The HXS 
light curve is plotted as a histogram to reveal its temporal 
correspondence with the HXR and NoRP data. Sometimes 
the HXS time bin could contain the emission integrated 
over two peaks seen by HXT and at 17 GHz (for example, 
this happened at about 04:03:47). Fortunately, the tempo- 
ral variations of HXR and NoRP emissions during the main 
peak resemble the HXS light curve. 



550 



S 500 



450 



400 




550 



550 



540 



£ 530 



520 



510 



600 650 700 

solar X [arcsec] 



750 



500 




560 



570 



-| — 1 — 1 — 1 — 1 — |- 

580 590 
solar X [arcsec] 



600 



610 



Fig. 3. a) Contours of 17 GHz (NoRH, 04:04) flare sources 
(black solid: Stokes I, 50%, 70% ,90% of the maximum, 
white sohd: Stokes V, RCP, 70%, 90% of the maximum, 
white dotted: Stokes V, LCP, 70%, 90% of the minimum) 
superimposed on an MDI magnetogram (04:48:01.61, light 
areas represent N-polarity, dark areas show S-polarity) . The 
axes show hereafter arc seconds from the solar disk centre, 
b) Enlarged fiare site denoted in panel (a) by the broken 
frame. Contours of 17 GHz (NoRH, 04:03:40) Contours lev- 
els are the same. Black thick contours correspond to the 
HXR source (Yohkoh/HXT/H; 04:03:40.25; 20%, 50%, 80% 
of the maximum). Thin solid gray contours mark N-polarity 
[levels 600, 800 G], and thin white contours mark S-polarity 
[levels -600, -500 G]. 



The time profiles of the HXR and 17 GHz emissions have 
a close temporal correlation (Fig. [5]). Nevertheless, there are 
some delays of HXR emission relative to microwaves. They 
amount to a fraction of se cond at high frequenci es and reach 
1-2 s at low frequencies (jAltvntsev et al.ll2008[ ). 
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2.2.2. MW and HXR images 

The flare was triggered by an interaction between two loops, 
one of which was a small, newly emerging loop, while the 
other was a large overlying loop (Chandra et al. 2006^ . The 
two loops formed a 'three-legged' configuration in which the 
magnetic field had a 'bipolar + remote unipolar' structure. 
The footpoints of the small loop (main flare source) can be 
seen in HXR/H and MW images (Fig. whose position 
is denoted by the white dotted frame in Fig. [3^. The mag- 
netogram taken at 04:48:30 was compensated for the differ- 
ential rotation to 04:03:40 (the most prominent peak time). 
The magnetic field strengths were projection-corrected us- 
ing the zradialize SolarSoftware routine (Fig. [3]). The length 
of the loop visible in an FIT 195 A image is 10000 km 
according to lAltvntsev et all (|2008[ ). During the flare, the 
NoRH beam size was about 20.6" x 17.3" at 17 GHz and 
9.5" X 3.9" at 35 GHz. 

The bulk of the MW emission at 17 GHz during the 
main peak was generated from the small loop. At the onset 
of the burst an LCP source located in the N-polarity mag- 
netic field (Fig.[3l3) was observed. An RCP source appeared 
north of the LCP source at 04:03:27 in the S-polarity mag- 
netic field, meaning that they both were polarized in the 
sense of the 0-mode. 

The magnitude of the field in the spot labelled "N" 
(solid thin gray contours in Fig. [3}d) reached 800 G and for 
the S-polarity field labelled "S" it reached 600 G. The large 
loop had footpoints in the main and remote (RS) sources 
(Fig. |3^). RS appeared at 04:03:51 and had a right circular 
polarization at 17 GHz. The right polarization in the main 
source disappeared after 04:03:52. 

At the onset of the flare (04:03:12) a single compact 
source was only seen in all energy channels of Yohkoh/HXT. 
HXR images obtained with a 4-s cadence in the L chan- 
nel showed that the source expanded and got elongated by 
04:04:05. In the H channel, the compact source transformed 
to a loop-like one by 04:03:25, and by 04:04:05 it became 
compact again. The distance between the centroids of the 
radio and HXR sources was about 10". 

For calculations of HXR spectra one needs areas of emit- 
ting sources. We estimated the sizes of the sources obtained 
in all four Yohkoh/HXT energy channels by fitting them 
with ellipses at half magnitude for the positions as if these 
sources would be located at the solar disk centre to exclude 
projection effects. The estimated sizes of the HXR sources 
are 16" x 12.6", 18" x 12.2", 15.4" x 12", and 12" x 16.8" 
for the L, Ml, M2, and H channel, respectively. 



There are the two basic expression to describe the prob- 
lem: a kinetic Fokker-Planck-Landau equation for beam 
electrons, and the Ampere law for the neutralization of elec- 
tric currents formed from the direct and returning beams. 
A distribution function / of beam ele ctrons is gover ned by 
Fokker-Planck, or Landau, equation ()Landaulll937D : 



— + Vcosa— 
ot az 



e£V cos I 



nirV dcos( 



df e£ s\v? a df 
"dE 
dl 

at 



dl 

dt 



(1) 



where z is a linear depth measured from the top of a coro- 
nal loop, E and a are the electron energy and pitch angle, 
respectively, V is the electron speed, and e and mo are 
the electron charge and mass, respectively. The first term 
on the left-hand side describes variations of the distribution 
function in time (i), the second term describes its variations 
with depth (z), the third and the fourth terms reflect energy 
losses owing to Ohmic heating and pitch-angular diffusion 
under the presence of a self- induced electric field £. The two 
terms on the right-hand side {df /dt)co\i and (9//i9i)magn 
describe the energy of a particle and the pitch-angle diffu- 
sion caused by scattering on ambient particles (collisional 
integral) and in a converging magnetic field, respectively. 

The colli sional integral is tak en in the linearized form 
suggested bv lDiakonov fc Somovl ([1988), considering varia- 
tions of collisional energy losses and anisotropic scattering 
with pitch angle diffusion for a given kinetic temperature 
of the ambient plasma and distribution function / of beam 
electrons as follows 



5/ 
dt 



dv 



fceTc df 
rrir dv 



ME) 



d 



dcos( 



sin^ a 



vf 

df 
9 cos a 



(2) 



Here fee is the Boltzmann gas constant, rc(^) is the 
ambient electron temperature at column depth ^ and i^(-E) 
is the frequency of collisions taken as follows 



lyiE) 



4fc V27r 



where n is the ambient plasma density, E is the beam 
electron energy, A(^) is the Coulomb logarithm for collisions 
of beam electrons with the ambient ones for a column depth 
^, and the parameter k is taken in the form (jEmslidl 19781) 



3. Model simulations 

3.1. Description of the model 

Let us consider the time-dependent Fokker-Planck ap- 
proach for precipitation of high-energy electrons injected 
into a cold hydrogen plasma conf ined in a convcrg- 
ing magnetic field structure (Leac h fc Petro sian 1981; 
ISiverskv fc Zharkovall2009l ) while producing a self-induced 
electric field, which forms a return current from the 
ambient plasma an d beam electrons (jKnight fc SturrocS 
Il977at lEmslid[T980f ). We take account of collisions, Ohmic 
losses and pitch-an gle anisotropy in scattering on ambi- 
ent particles (Diakon ov fc Somovil98l lMcClementill992al: 
IZharkova et~al.„201Ql . 



k^2x + [l- x) 



A" - 
A 



where A" and A' are the Coloumb logarithms for collisions 
of beam electrons with the ambient ions and neutrals, re- 
spectively. 

The linear coordinate z is replaced by the col- 

z 

umn density^(2;) = / nis) ds ()Diakonov fc Somovl 119881 : 



ISiverskv fc Zharkovall2009l) . The effect of the magnetic field 
convergence variations with c olumn density is described by 
the semi-empirical expression ISiverskv fc Zharkoval (j2009t) 
with the characteristic column depth of 10^° cm~^, allow- 
ing to match the observations of magnetic field magni- 
tudes at the photosphere and chromosphere (jKontar et al.l 
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2008fl. This semi-empirical form was furtlier modified by 



Siverskv &: Zharkoval (|2009[ ) for the magnetic field distr 



bution with a linear depth, required for the simulation of 
MW emission. 



3.1.1. Self-induced electric field 

The problem of electrostatic electric field carried by beam 
electrons and the return current formed by the ambient 
plasma was discusse d in the pa st decades by many authors 
(see, for example Knight fc _S turrock 1977b; E mslic 1983, 
and the references therein) . IBtowu fc Binghaml (j 19841 ) con- 
cluded that a return current can be established electro- 
statically because of large radii of electron beams inside flar- 
ing atmospheres that implies very long resistive timescales 
along the finite beam lengths leading to negligible inductive 
effects compared to the electrostatic ones even for plasma 
with anomalous resi stivity. The fu l l elec tro-magnetic ap- 
proach carried out bv lvan den OordI (119901 ) to study of elec- 
tric and magnetic fields induced by precipitating electron 
beams and their interaction with electrostatic field carried 
by these electrons allowed to conclude that the both elec- 
tric fields: electrostatic and induced ones co-exist in a flar- 
ing atmosphere and are neutralised by the electrostatic or 
solenoidal response of the ambient plasma. 

These t wo plasma responses ar e found to act inde- 
pendently (jBrown fc Binghaml 11984 Ivan den Oo^ |1990D 
forming from one side a return current from the ambi- 
ent electrons moving in the opposite direction to pre- 



trie field (iDiakonov fc Somov 


19881 Zharkova et al. 


19951: 


IZharkova fc Gordovskvv 12006 


; iBattaglia fc Bend 


20081). 



and from the other side by ambient electrons compen- 
sating for solenoidal electric (and magnetic) fields which 
can be oft en observed as longitudinal transient mag- 
netic fields (iKosovichev fc Zharkoval l200il : IZharkova et al.l 



l2005t ISudol fc HarvevI 120051 ). Therefore, for the p urpos e 



of this study, similarly to Diakonov fc Somovl ()1988D ; 



^ riy , ^ 

IZharkova fc Gordovskvvl (|2006l ). we consider only electro- 
static part of the electric field induced by beam electrons 
and assume that the electric current carried by beam elec- 
trons in the ambient plasma serving as the conducting me- 
dia is compensated by the return current formed from the 
ambient electrons. This allows us to estimate a magnitude 
of the electro-static electric field indu ced by the electron 
beam in the ambient plasm a as follows (jPiakonov fc Somovl 



The electric field induced by beam electrons carries a 
dual role: firstly, it induces (within the timescale of the 
double collisiona l time) the electro-m agnetic field in the 
ambient plasma (|van den OordI [l990l ) and sets up the re- 
turn current in the ambient plasma, and secondly, turns 
the precipitating beam electrons to the direction opposite 
to their original one, returning them back to the source 
and, thus, reducing the magnitude of the return current 
from the ambient elec trons by the factor of two or higher 
(jZharkova et al.ll2010D . 



3.1.2. The initial and boundary conditions 

It is assumed that at the moment of injection there are no 
beam electrons in the atmosphere, i.e. /(£^,/i, s, 0) — 0. 
The distribution function of beam electrons on the top 
boundary ^ = ^min is assumed to have power law distri- 
bution in energy in a range from Smin to i?max with a spec- 
tral index 7 -I- 0.5 (that corresponds to the particle energy 
flux varying with energy as a power-law with the index 7 
(Syrovats kii fc Shnieleva 1972)) and normal distribution in 
the cosines /i of pitch-angle a {p, — cos a) with the half- 
width dispersion A/i as defined below: 



fiE,fi,t)\ 



A^2 



Uit), 



beain m tne ambient piasni a 
[19881 IZharkova et al.l[l995l) : 



for E^in <E< E„,,^, fi>0, t>0, 
0, elsewhere. 

. . 

This condition represents an electron beam injected 
with a power-law dependence in energy and a narrow dis- 
persion (A^ <^ 1) in the pitch angle centered at a = (or 
H = 1). U(t) is a temporal profile of the electron beam injec- 
tion, which denotes the initial beam flux variations during 
a required time interval, accepted to be equal to unity (a 
steady injection). 

The coefficient A is found from the normalization con- 
dition for an electron distribution function as follows: 

^ - ' 

/ E^/^dE J fiE,^i,Uin)dii 

where Fq is the initial energy flux of accelerated electrons 
on the top boundary, normally derived from observations 
and set up as a free parameter in simulations. 



£{0 



^(0 ^(0 



v{E)dE / fi^,E,^i)^id^i, (3) 



where cr(^) is the classical plasma conductivity at a given 
temperature of the ambient plasma defined as (Spicer, 
1977): 



1 7.28 X IQ-^X , / 3 



7^1.5 

c 

7.6 X 10-1* 



In 



V2e3 



'^_B±c_ 

nn 



{1-x: 



X 



T 



0.5 



(4) 



where X is the ionization degree of the ambient plasma. 



3.1.3. Method of solution and accepted parameters 

The set of equations ^ and ([3]) defines the electron beam 
distributions in a flaring atmosphere at every instant of pre- 
cipitation. For their solution we use the summary approx - 
imation method described by ISiverskv fc Zharkoval ()2009() 
splitting the time interval into three bits and solving si- 
multaneously the temporal equations for electron distri- 
bution changes in depth, energy and pitch angle cosines. 
The solutions are sought for a electrons with relativistic 
correction precipitating in a quasi-stationary regime, e.g. 
when the electron distributions in depth and pitch angles 
do not change in time. This is normally achieved within 
0.07-0.2 s after the beam onset when the flows or pre- 
cipitating and returning e lectrons form a steady circuit 
(ISiverskv fc Zharkovall2009l) . 
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4. Results of simulations 

4.1. Estimation of beam parameters 



cm" 



10^ 



10" r 



10' 



f , cm 

Fig. 4. Height profiles of the plasma parameters (electron 
density and temperature vs. column density) which are used 
in the simulations. 



The simulations are carried out for electrons in the en- 
ergy range from 12 keV to 1.2 MeV. Electron distribution 
functions at higher energies (up to 10 MeV) were obtained 
by extrapolation of the mentioned simulations by using 
power- law fit. Comparison with the results of simulations 
covering a wider energy range has shown that such an inter- 
polation provides very high accuracy, while greatly reducing 
computation time for the relativistic FP equation. 

The following precipitation models are used in our sim- 
ulations: model C including only Coulomb collisions with 
the charged particles and neutrals of the ambient plasma; 
model C-l-E including the collisions and Ohmic losses in a 
self-induced electric field; model C-l-B including the colli- 
sions and scattering in converging magnetic field and model 
C-l-E-l-B considering all the above factors (collisions, elec- 
tric field and magnetic field convergence). 

In a flaring loop, both the electron density and tem- 
perature depend on height. These dependencies affect the 
plasma conductivity (jH) and coUisional integral ([2]) and 
thus must be taken into account. In this work, we use the 
temperature and den sity profiles derived from t he hydro- 
dynamic simulations (jZharkova fc Zharkovll200"7t ) shown in 
Fig. m We assume that the boundary of the transition re- 
gion is located at the depth where the column density f 
equals 10^" cm~^, and the profiles demonstrate a sharp 
change of the plasma parameters around this depth. 



lAltvntsev et al.l (|2008l) evaluated some plasma parameters 
from observations, so that the peak injection rate of elec- 
trons above 10 keV was found to be 8.5 x 10'^^ electrons s~^ , 
the total number of emitting electrons in the radio source 
was 3.83 X 10"^^ electrons, the background plasma density 
was 3 X lO^^cm^'^, the single power-law electron spectrum 
index 7 was 2.4, angular scale of the MW source was 6", the 
magnetic field value at the photosphere near the footpoints 
of the flare loop was 170-340 G. 

The authors (jAltvntsev et al.l[2008[) received the best fit 
for the following beam parameters: the characteristic life- 
time of the emitting electrons in the radio source ti=0.45 
s, viewing angle between the line of sight and the direction 
of the magnetic field at the source 0=80°, a beam precipi- 
tating along the direction of /io = 0.5 with the pitch angle 
dispersion at the half-width corresponding to A/i — 0.35. 

The current estimations of the HXR emission observed 
by Yohkoh (see section 12. 2p gave the total energy flux- 
about 1.5 X 10^^ erg/s/cm^ for energies from 10 keV to 
100 keV a nd 7 = 2.4 that was v ery close to the one es- 
timated by lAltvntsev et all ()2008[) . However, we extended 
the HXR photon energy range to 200 keV and for this en- 
ergy range (80-200 keV) we derived the spectral index to 
be 7 = 3.12. Hence, we have at hands the double power law 
photon spectrum with a break energy about 80 k eV. 

As shown by Zharkova fc Gordovskvvl ()2006[ ): 
ISiverskv fc Zharkoval (|2009t ) such double law HXR 



photon energy spectra are formed by combined collisional 
and Ohmic energy losses during precipitation of a single 
power law electron beam with a spectral index of beam 
electrons equal to that of the upper photon energy range 
(80-200 keV). Ohmic losses lead to flattening of he HXR 
photon spectrum at lower energies, which is higher for 
beams with higher initial energy fluxes and spectral 
indices. ^^^^^^^^^^^^^^^^ 

Also it was sh own (I Zharkova fc Gordovskvvl 120061 : 
IZharkova et al.l[201lt ) that, in order to calculate a total en- 
ergy flux of beam electrons accountable for the observed 
HXR photons, one needs to prolong the HXR photon spec- 
trum at higher energy (with index 3) to a lower energy 
range, to account for a single power law electron beam caus- 
ing this emission. Since the higher energy spectral index 
is higher than the one calculated with the spectral index 
(2.3)derived at the lower energy range we obtain the new, 
enhanced, total energy flux F ~ 10^^ erg/s/cm^ which ac- 
counts for both collisional and Ohmic energy losses by beam 
electrons. 

This difference in spectral indices and energy fluxes 
for lower and higher energy bands resembles very closely 
the HXR photon spectra produced in the models of elec- 
tron beam precipitation with the self-induced electric 
fleld which causes their flatten ing towards lower energies 
(jZharkova fc Gordovskvvl l2006l ). This indicates a need to 
consider the model of electron beam precipitation taking 
into account Ohmic losses, in addition to collisions, and for 
the beam parameters to adopt those derived for the higher 
energy range. 

In addition, because of the announcement in 2008 by 
MDI team about the problems with the magnetic field cali- 
bration leading to twice or smaller magnitudes of the mea- 
sured magnetic field with the previous one, we needed to 
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correct the measured magnetic field in the footpoints. With 
new calibration the magnetic field magnitudes at the pho- 
tosphere are increased to 600-800 G compared to the mag- 
nitudes of 170- 3 40 G accepted in the previous study by 
lAltvntsev et al.l ()2008f ) [the importance of the correction 
of the magnetic field strength to rec oncile MW and HX R 
spectra was recently demonstrated bv lKundu et alJ (|2009[ )] . 
Also, the calculations of HXR and MW emission in the cur- 
rent paper are carried out for the the visible HXR source 
area of about 150 arcsec^ (see section l^?^ . The fact that 
in this flare there is a presence of broadband MW emission 
in a wide frequency range (1-80 GHz) simultaneously with 
HXR emission from 10 to 200 keV allows to assume that 
they both are produced by the same populations of elec- 
trons, one needs to accept the energy for beam electrons 
ranging from 12 keV to 10 MeV. 

Another improvement of the previous simulation model 
was related to consideration of a flaring atmosphere highly 
inhomogeneous in atmospheric depth that was contrary 
to the homogeneous models of both ambient and beam 
ele ctrons considered in each of the two footpoint sources 
by lAltvntsev et al.l ([2008). The density and temperature 
variations are derived as a result of hydrodynamic re- 
sponse to the heating caused by injection of beam electrons 
(jZharkova fc Zharkovll2007[ ). This model produces signifl- 
cant depth variations not only in the physical conditions in 
the ambient plasma but also of the beam electrons them- 
selves if different energy loss and pitch angle scattering 
mechanisms are considered. 



4.2. Simulated HXR and MW emission 

Both HXR and MW emissions are assumed to be pro- 
duced by the same population of power law beam elec- 
trons with a spectral index of 3 ranging energies from 12 
keV to 10 MeV as derived in section 14.11 These beam 
electrons are steadily injected into a flaring atmosphere 
along the pitch angle zero {cosfi — 1) with a normal dis- 
tribution are getting scattered by the ambient particles, 
self-induced electric and converging magnetic fleld. The 
electron distribution functions for different models of en- 
ergy l osses (C, CE, CB and CEB) are discussed in de - 
tafl in'Siverskv fc Zharkova' ('2 009t ) : IZharkova et all (|2010D : 
iKuz nctsov & Zharkova (2010). This defines the electron 
beam dynamics and pitch-angle distributions at every pre- 
cipitation depth. 

The main issue in electron distributions is a formation 
of returning electrons either by magnetic field convergence 
or by self-induced electric field of precipitating beam elec- 
trons. These two beams, precipitating and returning ones, 
are producing either HXR emission in scattering on the am- 
bient plasma particles or slowing down in the electric field 
and MW gyro-synchrotron emission from their gyration in a 
strong converging magnetic field. Electrons emit HXR pho- 
tons in the directions downwards (to the photosphere with 
pitch angle a — — 90°) and upwards (to the observer at 
the Earth, with a = 90 — 180°) if the loop stands near the 
solar disk centre. If the loop is located outside the central 
zone, the viewing angles for HXR and MW observations 
will be amended by the angle of the loop position (latitude 
an d longitude) on the so lar disk which described in detail 
bv IZharkova et all (|2010[ ). 



4.2.1. HXR emission and directivity 

We calculate photon spectra of HXR emission integrated 
over all atmospheric depth, polarization and directivity nor- 
malized on the average intensity over all angles for the rel- 
ativistic cross-sections with pitch angle and viewing angle 
dependence as described below. 



A 



27^7^2 



X OO 1 

dC J v{E) AE J /(e,^;,^) d^^ 



hf 



2tt 

I 



{hiy,e,E,fi,ip) dip, 



(7) 



where A is the cross-section area of a magnetic tube, R 
is the astronomical unit, / and Q are the Stokes parame- 
ters related to the intensity and linear polarization of the 
emission. 

Then a degree of linear polarization can be defined as 
follows 



(8) 



and the distribution of any emission in the angles 9 can be 
described by a directivity 



D{9) 



m 

(!) ' 



(9) 



where (/) is the emission intensity averaged over all the 
angles. The three-dimensional integrals over the electron 
velocity in ([7]) are calculated with Monte-Carlo method. 

The simulations of HXR photon spectra for different 
precipitations models (C, GB, CE and CEB) including 
those emitted downwards (to the photosphere) and upwards 
(to the observer) and for different magnitudes of viewing 
angles are plotted in Fig. [5] from the top to the bottom, 
respectively. On top of the simulated energy spectra we 
over-plotted by a solid line the energy spectrum deduced 
from the observations (see section [2?2t . 

The plots of HXR emission simulated for different direc- 
tions of electron propagation (downwards with pitch angle 
cosines fj, > and upwards with /i < 0) reveals that in the 
corona the majority of electrons moves in the downward di- 
rection while only a twice smaller number moves upwards. 
In the chromosphere, it is to the contrary: most electrons 
move upwards and only a smaller fraction of them keeps 
moving downwards. Thus, for this flare both effects: albedo 
and Ohmic losses have to be very significant at upper (coro- 
nal) precipitation depths. 

The spectra calculated for GB model (collisions and 
converging magnetic field) sho w the highest intensity com- 
pared to CE and CEB models ijZharkova et al.ll2010l ). This 
is because the electric field induced by precipitating elec- 
tron beam prevents the particles from reaching deeper lay- 
ers where the bulk of electrons (with lower cutoff energy) 
can emit HXR photons, thus reducing the intensity. The 
CE model provides a smaller intensity than for model GB 
while the intensities for CEB model are much smaller than 
the either of the above for CB or CE. This happens be- 
cause of the combined effect of the two factors (convergence 
and electric fields) which significantly reduces a number of 
particles being able to precipitate into deeper atmospheric 
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Fig. 5. Intensity (in arbitrary units) of the HXR emission 
spectra calculated from equation ([7]) for different energies, 
models and propagation directions (0 to 90° for downward 
emission and from 90° to 180° for upward emission). Dash- 
dotted line: emission from upward propagating particles; 
dotted line: from downward propagating particles, dashed 
line: total emission (downward + upward); solid line: ob- 
servational data multiplied by the corresponding factors to 
match the simulated absolute intensities. 



layer s where the bu l k of electrons can emit HXR pho- 
tons. IZharkova et all (|2010D have shown that for F = 10^° 
erg/s/cm^, the results were opposite: the CE model pro- 
vides a higher HXR intens ity than the CB mo del. This 
confirms the conclusion by (jZharkova et al.ir2010l) that for 
F = 10^^ erg/s/cm^ the effect of a self- induced electric field 
is much stronger than for a weaker beam. 

It can be also noted that the observed HXR photon 
spectrum has the best fit by the HXR emission emitted 
upwards plus downwards (the full albedo effect) under a 
viewing angle of 180° by an electron beam with given pa- 
rameters for the CE (collisions and electric field) and CEB 
(collisions, electric field and converging magnetic field with 
convergence equal to three) precipitation models. The HXR 
spectrum observed in this fiare reveals a noticeable flatten- 
ing towards lower energies which is better reproduced by 
CE model indicating a significant effect of the self-induced 



electric field in the energy losses by beam electrons. The 
observed and simulated energy spectra reveal the best fit 
for the viewing angle of 180° (the direction towards the 
observer looking from the top). Given the flare location 
at the location (N27W42) one can assume that the flaring 
loop emitting HXR photons is tilted towards the observer 
by 10° - 40°. 

At the same time the directivity of HXR emission is 
found to range between 2 and 3 for the electrons propa- 
gating downwards and about 0.5 for the electrons moving 
upwards. This points out to clear anisotropy of the electron 
beam producing this HXR emission with domination of the 
downward moving electrons. In order to see this emission 
from the top where the observer is placed in our models, 
HXR photons have to be reflected by the photosphere, e.g. 
the photosph e ric alb edo effect should play a significant role 
iKontar et al.l ([200^. In the present study we have assumed 
100% albedo coefficient for the emission emitted downwards 
while in reality this can have more complic ated dependence 
on energy of electrons (jKontar et al.ll2006h and their pitch- 
angle distribution. However, this effect requires to consider 
strong electron anisotropy which we described above, con- 
trary to the i s otrop ic electron distributions considered by 
iKontar et al.l (|2006[ ) that will be considered in the forth- 
coming paper. 

4.2.2. MW emission 

We now calculate the parameters of MW emission. The 
model of the emission sourc e is similar to that used by 
iKuznetsov &: Zharkoval (j2010[ ) : we assume that the param- 
eters of the coronal magnetic tube (such as the plasma den- 
sity, magnetic field, and the parameters of the accelerated 
electrons) depend only on the coordinate z, a linear dis- 
tance along the tube. Also, we used a stratification of MW 
emission across the layers in vertical direction (correspond- 
ing to different heights in the atmosphere). These layers 
are assumed to be quasi-homogeneous sources, so that the 
emissivity and absorption coefficients within each layer are 
assumed to be constant. As result, the total MW emission 
emitted from the whole atmosphere is equal to a sum of the 
contributions by all layers written as follows: 

2m ax 

" " 7^2 J x,{z) [ \ ^ ^ 



Here /o- is the emission intensity (observed at the Earth) 
of the magneto-ionic mode a, D and L are the visible diam- 
eter of the magnetic tube and the source depth along the 
line-of-sight, respectively, and TZ is the astronomical unit. 
The equations for the gyrosynchrotron plasma e missivity 
jg- and abs orption c o efficie nt Ko- are given by e.g. iMelrosd 
(|1968D and iRamatvl (|1969D . The polarization degree is de- 
fined as 

-^x ~ Iq , . 

V= r ,r , (11) 

where Iq and Ix are the intensities of the ordinary and 
extraordinary modes, respectively. The directivity of MW 
emission is calculated for different propagation (or viewing) 
angles) , similar to that of HXR emission (see equation ^ 
with the intensity calculated by using equation (jlOp . 
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Fig. 6. Directivity of the HXR emission calculated from equation ([9]) for electron energies of 30 keV, 50 keV, 100 keV 
and 200 keV for different models of electron energy losses and propagation directions. Solid line for CE model, dotted 
line for CB and dashed line for CEB model. Propagation direction indicates the viewing angle for the observer looking 
from the top and varies from to 90° for downward emission and from 90° to 180° for upward emission. 



Unlike the previous paper (jKuznetsov &: Zharkoval 

l2010t) . in the present study we use the temperature and 
density of a flar ing model der ived from the hydrodynamic 
simulations ( Zh arkova fc Zhar kov 2007) with the distance 
z from the injection point varying from to ^max — 10 000 
km, the thermal plasma density varying from 2 x 10® cm^"^ 
to 2 X 10^'^ cm^'^ (see Fig. H]), and the magnetic field 
strength oi B ~ 780 G at the characteristic column depth 
of = lO^'^ cm^-^ (or z — Zmax)- The parameter Zmax corre- 
sponds to the boundary of the transition region, since only 
the coronal part of the loop makes a contribution into the 
MW emission (in the deeper layers, the plasma density is 
too high). 

In the models with converging magnetic field, the field 
strength varies with depth from the top of the emission 
sourc e to this characteristic depth (jSiverskv fc Zharkoval 
l2009 lh , which is another factor (in addition to the variations 
of a distribution function) affecting the MW emission pa- 
rameters; we used the models with the convergence factors 
-Bfootpoint/^^top = 2 and 3 which provided the very reason- 
able fit to the observations (see sections [4.3.11 and 14 . 3 . 2p . In 
the CE model, variation of the magnetic field with height 
(with the convergence factor 3) was considered when cal- 
culating the plasma emissivity and absorption coefficient, 
but the effect of converging magnetic field on the electron 



distribution was neglected. Also, we assume that the loop 
width is 13 = 10 000 km and the source depth along line- 
of-sight is L = 10 000 km; these parameters agree with the 
imaging observations and provide the best agreement of the 
calculated MW spectra with the observed ones (see below). 

MW emission was simulated with radiative transfer ap- 
proach for vertically stratified layers for the following beam 
precipitation models described in section |3] (CB, CE and 
CEB) and for different viewing angles (see Fig. [7]). It can 
be noted that for the viewing angle 6 — 110°, the main fac- 
tor affecting MW emission is a magnetic field convergence. 
MW intensities have a rather fiat maximum at about 20 
GHz approaching the magnitude above 10"^ sfu for the mod- 
els with a magnetic field convergence (CB and CEB mod- 
els) and reducing the maximum below 10^ sfu at 10 GHz 
without it (CE model). In the models including a magnetic 
field convergence without (CB) or with electric field (CEB), 
the plots of MW intensity and polarization are rather close 
with the intensities decreasing for a convergence factor of 3 
compared to 2. The models with the converging magnetic 
field (CB and CEB) provide much higher emission intensity 
than the model without this factor (CE). 

The effect of a self-induced electric field is relatively 
weak and is visible only if the magnetic field convergence 
is equal to zero, e.g. the model with the return current 
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Fig. 7. Intensity (equation (fTO|) ') in solar flux units (sfu) (upper plots) and degree of polarization (equation (|TT|) ) (lower 
plots) of the MW emission spectra calculated for different precipitation models and propagation directions: for the viewing 
angles of 110° (left column plots); 140° (middle column plots) and 170° (right column plots). Numbers in the model 
abbreviations define the factor of magnetic field convergence, e.g CB2 or CEB3 while the number in brackets for CE(3) 
denotes simulations for DF calculated for CE model with magnetic field added in the MW intensity calculations. Note 
that propagation direction indicates the viewing angle for the observer looking from the top and varies from to 90° for 
downward emission and from 90° to 180° for upward emission. 



(CE) provides a lower peak intensity and a steeper inten- 
sity decrease towards higher frequencies than the collisional 
model with magnetic convergence (CB). The MW intensity 
at higher energies are power laws with spectral indices of 
about 2.4 for CEB models and higher than 4 for CE model. 
The MW frequency distributions towards lower frequen- 
cies reveals also power laws with lower spectral index of 2 
revealing strong harmonic structure from frequencies of 8 
GHz towards zero. 

The polarization produced by beams for this viewing 
angle of 110° in CE models has always negative sign and 
ranging in a few percent interval while the polarization pro- 
duced by CB or CEB models with convergence factor of 2 
or 3 is positive (X-mode dominates) for higher frequen- 
cies approaching a few percent and crossing zero between 
10 and 11 GHz after which it becomes negative (0-mode 
dominates). The harmonic structure in MW polarization is 
more pronounced than that in the intensity, and it increases 
with a growth of magnetic field convergence. 

For a viewing angle 9 — 140°, the effect of a self- induced 
electric field becomes much more significant than those of 
the converging magnetic field: the models with collisions 
and self-induced electric field (CEB) provide higher MW 
intensity than the models with collisions and magnetic field 
convergence (CB). The maximum MW intensity increases 
for the models with electric field and magnetic field conver- 
gence (CEB) compared to those at 110° and frequency of 
maximum shifts from 11 to 12 GHz for 140°. The power law 



spectrum at high frequencies has spectral indices ranging 
from 2.4 to 2.8 for CEB3 and CEB2 models, respectively. 
While the lower energy part of MW intensity resembles the 
power law distributions found for 110° with wave- like os- 
cillations with slightly larger amplitudes. 

The self-induced electric field has a significant effect on 
the MW polarization which now changes the sign for all 
models, and the frequency of this change shifts to the mag- 
nitudes below 10 GHz. The MW polarization degree for CB 
models is higher than those with electric field (CEB) ap- 
proaching 50% at 11 GHz for CB2 model with convergence 
factor 2 and reducing to 40% for CB3 with the convergence 
3. The polarization for models CEB becomes slightly lower 
that for CB with rather flat distribution towards higher 
frequencies and strong wave-like oscillation towards lower 
frequencies with zero-points between 9 and 10 GHz. 

For a viewing angle 170° the effect of the self-induced 
electric field is still noticeable for all the models, although 
the models with collisions, converging magnetic field and 
self-induced electric field (CEB) provides slightly higher 
MW intensity than the models with collisions and electric 
field (CE) only.However, the maximum MW intensity is sig- 
nificantly (by order of magnitude) reduced for all the mod- 
els with maxima ranging 7-10 GHz for CE and CEB mod- 
els, respectively. The power law spectra at high frequencies 
become much softer having the spectral indices 4-5. While 
the lower energy part of MW intensity resembles the power 
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law distributions found for 110° with harmonic structure 
less pronounced than for other viewing angles. 

The polarization also changes significantly compared to 
the other viewing angles approaching about 100% from fre- 
quencies 3-5 GHz for CBS and CB2 models, shifting to 4-6 
for the models CEB3 and CEB2, respectively. The polar- 
ization sign change is shifted to the frequencies of 2 GHz for 
CBS and 4 GHz for CEB2. The strong harmonic structure 
with very large amplitudes is still present at lower frequen- 
cies. 

These distributions of MW intensity and polarization 
can be understood after considering the combined effect of 
a self-induced electric field and converging magnetic field 
which turns the electrons moving upwards to such the pitch- 
angles that their distribution maximum occurs at about 
IX = —0.8 (Zharkova fc G ordovskyy .2006i) . Hence, the di- 
rection where the particles emit the MW (and HXR) ra- 
diation corresponds to the intermediate viewing directions 
(9 ~ 120° — 150°). As a result, taking into account a self- 
induced electric field (in the models GEB) results in the 
increase of a maximal MW intensity by a factor 3 (in com- 
parison with the model CB) and in a noticeable shift of the 
spectral peak from 10 to 20-30 GHz. 

The spectral index at high frequencies (for 140°) is -2.4 
for the CEB3 model, and -2.8 for the CEB2 model. This 
is close to the value of -3 derived from the observed MW 
intensity distribution in frequency. At low frequencies, the 
spectral index is about 3 that is higher than the observed 
one. This can mean that, though the source is large and 
inhomogeneous, in the numerical model the largest contri- 
bution of the MW emission comes from a compact region 
near the loop footpoint. 

4.2.3. MW emission directivity 

Now let us consider the MW directivity plots in Fig. [5] cal- 
culated for the same parameters as MW intensities in Fig.[71 

It can be noted that the directivity is also strongly de- 
pendent on the electron precipitation model applied in cal- 
culations. In all models there are two preferential direc- 
tions for emitting the MW emission: downwards and up- 
wards with the maxima position varying for different mod- 
els. For the precipitation models including only electric field 
(CE(3)) calculated for 20 GHz (left plot) and 30 GHz (right 
plot) the most of MW emission is emitted in the direction 
towards the photosphere with a maximum at about 30°. 
This maximum is twice higher than the maximum MW 
emission occurred in the direction of 135° (e.g. towards the 
observer from the top). 

For the models including collisions and magnetic field 
convergence the directivity of emission for MW emission at 
20 GHz (top left plot) still has the maxima at the same 
viewing angle, but the magnitudes of the maxima are re- 
duced in both directions, although, keeping the ratio be- 
tween the downward and upward emission at about 2 that 
is close to that for CE model. For the MW emission at 
30 GHz (top right plot) the ratio between downward and 
upward emission increases to 2.7 for CB2 model approach- 
ing 1.4 for CEB2 model. However, CB3 models show some 
isotropisation of MW emission between 30° and 130° where 
there is no preferential direction is found while inclusion of 
the electric field in CEB3 model reveals again some prefer- 
ential direction with a small ratio. 



Comparison of the CEB model simulations for differ- 
ent convergence factors (CEB3 in the bottom left plot and 
CEB2 in the bottom right plot) reveals that directivity 
is strongest for the model with convergence factor of 2 
with the downward-to-upward emission ratio approaching 
4, while for the convergence 3 it is only reach 1.5-1.8 for 
the lowest frequency. There is also a visible reduction of 
the directivity with the increase of frequency for any con- 
vergence although, the differences are much higher for the 
convergence factor of 3 (the left plot). These properties of 
MW directivity, in addition to intensities and polarization, 
help us to diagnose electron beam precipitation scenario in 
this flare. 

4.3. Fit to observations and general discussion 

The simulations of kinetics of a single electron beam with 
wide energy range from 12 keV to 10 MeV precipitating into 
a flaring atmosphere being heated by this beam via hydro- 
dynamic response is applied to simultaneously interpret the 
observed MW and HXR emission. 

4.3.1. Hard X-ray emission 

A comparison of the distributions versus energy of the sim- 
ulated and observed HXR photon emission (downward, up- 
wards and their sum) derived from YOHKOH instrument 
is presented in Fig. [5] for the two viewing angles of 90° 
and 180° and two sets of models: collisions and electric 
held (CE) and collisions, electric field and magnetic field 
convergence of 3 (CEB). It is obvious that the observed 
HXR photon spectrum cannot be fit by any models for 
the viewing ang l e of 8 0° derived from MW observation by 
lAltvntsev et al.l (|2008D . However, we can assume that HXR 
and MW emission come from different parts of the loop and 
this discrepancy between the derived viewing angles can be 
related to a loop curvature. 

Although the agreement between HXR observations by 
the YOHKOH payload and the proposed models (shown in 
Figure [S] in linear scale to amplify the differences) is much 
better for the simulations carried out for a viewing angle of 
180° (e.g. upwards along the magnetic field) for the models 
including electric field (CE and CEB). It is evident from the 
residuals that the fit is better for CE model (collisions and 
electric field losses) combining the upward and downward 
emission (the full albedo effect). 

This fit can indicate that, at first, albedo effects are 
rather important in the interpretation of HXR emission 
produced by well colli mated beams, in ad dition to isotropic 
ones proposed earlier (jKontar et al.l 12006) : and at second, 
the observed HXR emission is likely to come from the part 
of a fiaring atmosphere not affected by a magnetic field con- 
vergence, e.g. from column densities larger then the char- 
acteristic column density discussed in section |21 

4.3.2. MW emission and polarization 

The distributions in frequency of MW emission and polar- 
isation observed for this flare in 9 frequencies shown by 
asterisks were compared with the MW emission simulated 
with a radiative transfer approach in vertically stratified 
layers for different viewing angles (180, 140 and 110°) and 
for different electron precipitation models including coUi- 
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Fig. 8. Upper plots: directivity of the integrated MW emission spectra (equation ^ with the intensity from equation 
(fTO)l ) versus different propagation directions (see note in Fig.[7l) simulated for the models CB, CE, and CEB (indicated on 
the plots) with convergence factors 2 and 3 for 20 GHz (left plot) and 30 GHz (right plot). Bottom plots: MW directivity 
versus propagation directions calculated for different frequencies (indicated on the plots) with CEB3 model (left plot) 
and CEB2 model (right plot). CE(3) indicates the MW emission calculated for the distribution functions with CE model 
but with the magnetic field magnitudes at relevant depths corresponding to CEB3. 



Table 1. The statistics and significance coefficients for correlation between the observed and simulated MW intensities. 



Points 






Kendal 


Spearman 




Kendal 


Spearman 


x' 


of CEB 


0.90° 


0.99'' 


CEB3 " 


CEB3 


CEB3 " 


CEB2 " 


CEB2 


CEB2 " 


9 


3.49 


1.65 


0.944 


0.944 


3.86 


0.889 


0.950 


2.56 


7 


2.20 


0.87 


0.867 


0.943 


2.21 


0.733 


0.829 


0.331 


6 


1.61 


0.55 


1.000 


1.000 


1.87 


0.800 


0.900 


0.128 


5 


1.06 


0.30 


1.000 


1.000 


1.49 


1.000 


1.000 


0.031 



Notes. Xo for 90% significantce level from tables bv lChernoff fc LehmannI (|1954f ): ^''^ Xo for 99% significantce level from tables 
by IChernoff fc LehmannI (| 19541 ): Kendall tau correlation coefficient; Spearman correlation coefficient; x^ derived from 
comparison of the observed and simulated MW intensities for magnetic convergence 2 (CEB2) and 3 (CEB3). 



sions, electric field and magnetic field convergence as pre- 
sented in Fig. [TUlfor a viewing angle of 140° revealing the 
most reasonable agreement. 

The observed MW spectrum was created by using the 
NoRP data plotted in Fig. [T^ and the MW polarization 
variations were taken from the Fig. [T}d. The MW emission 



simulated for optically thin plasma used for fitting the ob- 
servations in the paper by lAltvntsev et al.l (|2008D is pre- 
sented on the plot by the dot-dashed line. The errors of 
MW emission measurements are shown by bars in Fig. 1101 
We calculated the standard deviation of the observed noise 
level (cr) for each frequency of MW observations by SBRS 



13 



Zharkova et al.: Diagnostics of the beam anisotropy 





4e7 
3e7 




\ model 

\ observations 




0.04 
0.00 


=_CE up J 




2e7 

1e7 





CE up ^ 




-0.04 






1b7 








0.04 


^CE down y ^ 




5e6 








0.00 











CE down^""""" — 




-0.04 




jnits 


4e7 






arbitrary unil 


0.04 


=_CE sum 


trary i 


2e7 




CE ,„m~'""--— — 


0.00 

-0.04 


r ~ 


itensity, arbi 


8e6 
6e6 
4e6 

2e6 





\ - 

CEB up^~^_____ i 


Residuals, 


0.04 
0.00 
-0.04 


LCEB up J 

j- -j 




2e6 










t!GE.^^jown ^ ^^^^^^^^^^ ^ y'' _i 




le6 








0.00 











CEB doJs~— 5 




-0.04 






le7 








0,04 


|_CEBsum J 




5e6 






CEB^^^^urn'^ ^ : 




0.00 
-0,04 



















50 100 150 200 250 300 ,00 

Energy, keV Energy, l<eV 



Fig. 9. Comparison of the simulated upward, downward and total HXR intensities with the observations (left plots) and 
their residuals (right plots) calculated for different beam precipitation models (CE and CEB) and the viewing angle 180° 
(along the magnetic field direction). 



and NoRP. The errors in MW intensity and polarisation 
were estimated then from a sum of the statistical (3cr) and 
instrumental errors. 

From a hydrodynamic model used to define the temper- 
ature and density variations with depth one can estimate 
the height of a semi-circular loop producing this flare to be 
about 10000 km. Then the distance between the MW and 
HXR intensity centroids can be estimated at about 10" that 
fits very well the locations of MW and HXR contours. 

The MW emission and polarisation simulated in sec- 
tion by taking into account radiative transfer effects 
reveal smoother parts at medium and higher frequencies 
and a harmonic part at lower frequencies. The smooth part 
at higher frequencies has a negative power law close to that 
observed while for lower energy part it has a positive power 
law with a rather steep slope which is much higher than the 
one observed. In the current paper we are mostly concerned 
with the smooth part of MW emission and polarization 
simulated at higher frequencies for convergence of 2 and 3 
(CEB2 and CEB3 and different viewing angles. While at 
lower frequencies we have accounted for Razin's effect with 
the consideration of depth variations of the ambient den- 
sity that improved the MW emission magnitude at lower 
frequencies. However, there are still some noticeable dis- 
agreements at very low frequencies between the observed 
and simulated MW emission as seen in Fig. \TU\ . These are 
possibly, related to the use of a simplified radiative transfer 
approach or some other effects which will require further 
investigation in the future. 

For a viewing angle of 11 0° , which in o u r mo del is 
equal to 70° in the model bv lAltvntsev et all (|2008[ ). the 
simulated MW intensity is not sensitive to the variations of 
the magnetic convergence factor and it is smaller than the 
observed one by more than 20% for both lower and higher 
frequencies. However, for a viewing angle of 140° the model 
simulations of MW emission and polarization are able to re- 



produce very closely the main features of the observed MW 
emission. These include the magnitude and the frequency of 
MW emission maximum close to the ones observed; the fre- 
quency of the MW spectrum maximum appearing at about 
17 GHz, a negative power law spectral index and a smooth 
decrease of MW emission towards lower frequencies with a 
spectral index about 2 (being still much higher than the 
observed one because of the simplified radiative transfer 
approach used and negative power law MW emission with 
the spectral index of 2.8 towards higher frequencies close 
to the observed one of 3). 

The MW polarization simulated for models CEB2 and 
CEB3 for the same viewing angle of 140° agrees with the 
observed polarisation at lower frequencies with CEB2 show- 
ing a closer agreement. Also the CEB2 model simulation of 
MW polarisation reproduces rather closely (9 GHz), within 
the limitations of our radiative transfer model, the observed 
frequency (7 GHz) of the polarisation reversal. While for 
higher frequencies the simulated polarisation is twice higher 
than the observations that can be the result of neglecting 
further scattering of MW emission on the ambient particles 
(jBastianlllQQSl lAltvntsev eraI1ll996D , which can smoothen 
significantly a polarisation degree of MW emission at higher 
frequencies. 

In order to quantify our fits in MW intensities, we car- 
ried out a few statistical tests for the likehood of each 
curves (CEB2 and observations, CEB3 and observations) 
by using the statistical SPSS package. The t ests include 
Kend all tau and Spearman co rrelation tests (|Isobe et al.l 
Il986l : iKendal fc Gobboiislll990t) appropriate for other than 
normal distributions which we have for MW emission. We 
also calculated coefficients ([Chernoff fc Lehmannlll953) 
for the following CEB2 and CEB3 sets: the full number of 9 
measurements, for reduced 7 measurements (excluding the 
two points at lowest frequencies) , then for 6 and 5 (reducing 
the measurements consequently by one from lower frequen- 
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cies), in order to avoid the discrepancies in MW emission 
at the lowest frequencies. The resuhs of these tests are pre- 
sented in Table 1. For comparison, we also include the crit- 
ical values of for given degrees of freedom corresponding 
to significance levels of 90% and 99%. 

It can be observed from Table 1 that the intensities 
for both models (CEB2 and CEB3) reveal strong positive 
correlation with MW observations for all 9 frequencies, but 
it increases to a full dependence (correlation coefficients 
approach 1.0) for the reduced number of measurements of 
6 and 5 for CEB3 or 5 for CEB2. Although, for 9 points 
the correlation is slightly better for the model CEB2. The 
calculated for each datasets (CEB2 and CEB3) from 
the full to reduced measurement numbers plotted in Table 
1 allow us to discriminate between the models. The model 
CEB3 produces of 3.86 for 9 points (or degree of freedom 
(DF) of 8) that is reduced to 1.49 for 5 points, or DF=4). 

The x^ coefficients for simulated and observed polari- 
sation data obtained for 8 observational points (excluding 
80 GHz), or 7 degrees of freedom, produced x^=l-33 for 
model CEB2 and 2.13 for model CEB3 that reveals the 
similar close (above 95% fit to observations for the model 
with convergence 2 because their x^ he between the critical 
values for the confidence lev els above 0.90 and below 0.99 
(|Chernoff fc Lehmand^l954^ . This analysis of both MW in- 
tensity and polarisation allows to conclude that model sim- 
ulations with any convergence fit reasonably well the MW 
observations up to the confidence level of 90%. However, 
the model CEB2 with convergence 2 fits the observations 
up to 95% for all 9 points and up to the confidence level of 
99% with a reduced number of 7 points and lower. 

Therefore, it is safe to conclude that with the electron 
beam parameters derived from HXR emission, the best 
agreement with the MW observations at higher frequencies 
is achieved for CEB models with a reasonably small con- 
vergence factor of 2, the magnetic field magnitude at the 
photosphere derived from MDI (B = 780 G) and the view- 
ing angle 140° (see Fig. ITUl) . The results of the simulations 
also allow us to estimate more accurately the inclination of 
the loop in the parts emitting the relevant types of emis- 
sion: HXR or MW. 

One can note that the observed MW emission from this 
flare better agrees with our simulated MW emission seen 
under the viewing angle of 140°, contrary to the viewing 
angle of 80° (equal to 100 ° in the current mod el) deduced 
in the earlier simulations (jAltvntsev et al.|[2008[ ) and to the 
angle of 180° obtained for HXR emission in the current 
paper. As we discussed above, the difference between the 
viewing angl es derived for MW em ission from our models 
and those bv lAltvntsev et all (|2008[ ) can be a result of the 
over-simplified model used in the latter and a limited num- 
ber of measurements. 

However, the difference between MW and HXR viewing 
angles found in the current paper can be explained within 
our model by the assumption that the flaring loop was not 
standing perpendicular to the local horizontal plane in the 
flare location but was slightly tilted (by about 40°) towards 
the solar disk center. By comparing the directivity of MW 
and HXR emission from Fig. [S] one can observe that the 
MW emission has the preferential direction of 30° down- 
wards and of 130-140° upwards. 

The downward emission exceeds by factor 2-4 the up- 
ward emission. If all these downward photons are fully re- 
flected by the photosphere (albedo effect) with the same 
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Fig. 10. Comparison of observed MW intensity in solar 
flux units (sfu)(top plot) and degree of polarization (bot- 
tom plot) with results of the simulations with CEB2 model 
for a magnetic field convergence of 2 (dotted line curves) 
and CEB3 model with conv ergence of 3 (solid line curves) 
and the model presented bv lAltvntsev et all (j2008|) (dash- 
dotted line) . In the both plots asterisks mark the observed 
spectra. 



properties (angles and energies) they can contribute to rea- 
sonable fit of MW emission for the viewing angle of 140° 
compared to 180° from HXR spectra, e.g. this difference in 
viewing angles can be just a result of the directivity effects 
on HXR and MW intensity from beam electrons. One can 
also speculate that the loop curvature in the observed flare 
could contribute to this angle difference for viewing MW 
and HXR emission to conflrm which one needs to extend 
our model to semi-circular loop that can be done in the 
future. 



5. Conclusions 

In this paper for the first an attempt is made to simulate 
HXR and MW emission and polarisation emitted from he 
flare of 10 March 2001 with the same population of elec- 
trons, in order to improve the fit of observations reported in 
previous studies. This is achieved by applying the Fokker- 
Planck kinetic approach to precipitation of electron beam 
with energy range from 12 keV to 10 MeV into a con- 
verging magnetic loop with anisotropic electron scattering 
on the ambient particles in Coloumb collisions and Ohmic 
losses. The theoretical HXR and MW emissions are then 
calculated by using these distribution functions for differ- 
ent factors of magnetic field convergence a nd viewing angles 
as described in our previous papers (Zhar kova et al.l 120101 : 
iKuznetsov fc Zharkovall201(l . 

The observed HXR photon spectra and frequency dis- 
tribution of MW emission and polarisation reveal the best 
fit for the FP models including the effect of electric field 
induced by beam electrons precipitating in a converging 
magnetic loop. Magnetic field strengths in the footpoints 
on the photosphere were updated with newly calibrated 
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SOHO/MDI data. The observed HXR energy spectrum is 
shown to be a double power law which was fit very closely 
by the photon spectrum simulated for the models including 
the self-induced electric field. 

The MW emission simulated for different models was 
compared with the observed distribution in frequency re- 
vealing that only the models combining collisions and elec- 
tric field effects with pitch angle anisotropy were able to 
reproduce closely the main features of the observed MW 
emission: peak magnitude and frequency, a negative spec- 
tral index of 2.8 versus 3 observed for higher frequencies, 
a smoother decrease of MW emission towards lower fre- 
quencies and the correct frequency of the MW polarisation 
reversal. There are still disagreements between simulated 
and MW emission observed at lower frequencies which are 
likely to be caused by simplified model of radiative trans- 
fer considered in the current approach and by not taking 
into account the emission from the ambient electrons in a 
vicinity of the flare. 

By estimating from the hydrodynamic model the height 
of a semi-circular loop to be about 10000 km and by taking 
into account the inclination of the loop in HXR emission, 
the distance between the MW and HXR intensity centroids 
was about 10" that agrees rather well the observations. 

HXR emission simulated for relativistic angle- 
dependent cross-section accounting for different directions 
of electron propagation (downwards and upwards) reveal 
that in the corona the majority of electrons moves in 
the downward direction with a twice smaller number 
moving upwards while in the chromosphere most electrons 
move upwards and only smaller fraction keeps moving 
downwards. Thus, for this flare the effects of electron's 
magnetic mirroring and Ohmic losses were significant. 

The observed HXR spectrum reveals a noticeable flat- 
tening towards lower energies below 50 keV indicating a 
significant effect of the self-induced electric field in beam 
electron energy losses. The observed HXR photon spec- 
trum was best fit by the model simulated for collisional plus 
Ohmic losses (CE) precipitation model of electron beam 
with the initial energy flux of 10^^ erg/s/cm^ and spectral 
index about 3 including the emission emitted upwards and 
downwards (the full albedo effect) in the direction of 180° 
(towards the magnetic field direction) that for this flare lo- 
cation indicates a loop tilt about 40° towards the solar disk 
centre. 

The observed MW emission distribution in frequency 
reveals a better flt for the model combining collisions and 
electric fleld effects with a moderate magnetic field conver- 
gence of 2 (CEB2). This electron precipitation model can 
reproduce closely within better than 90% confidence level 
for all 9 points of observations and up to 95% confidence 
level for 6 points excluding the three ones at lowest frequen- 
cies) the main features in the observed MW emission: the 
maximum magnitude and the frequency at about 17 GHz, 
a smooth decrease of MW emission towards lower frequen- 
cies with a spectral index about 2, the spectral index of 2.8 
for higher energy part of MW emission and the frequency 
of reversal of MW polarisation. 

The both models CEB2 and CEB3 show a reasonable 
fit to the observed polarisation magnitudes; although the 
model CEB2 shows a better fit to the frequency of the MW 
polarisation reversal close to the observed one within the 
90% confidence level for all 9 points of observations and 



up to 95% confidence level for 6 points excluding the three 
observations at lowest frequencies. 

In this study the closest fit for observed MW emission 
was achieved for the simulations with a weakly (factor 2) 
converging magnetic field seen at the viewing angle of 140°, 
contrary to 80° ( corresponding to 100° in our model) de- 
rived previously (jAltvntsev et all l2008\ While for HXR 
emission the closest agreement between observed and sim- 
ulated photon spectra was obtained for the model with col- 
lisions and electric field only with the emission seen at 180° 
(from the loop top). 

Within the limitations of our radiative transfer model, 
this difference is likely to refiect, at first, the inclination of 
about 40° of the fiaring loop towards the solar disk centre 
that accounts for the HXR viewing angle. At second, this 
difference can also indicate a difference in the preferential 
directivity of HXR and MW emission caused by electron 
scattering effects in a presence of self-induced electric and 
converging magnetic fields. In addition, a fiaring loop cur- 
vature in the upper atmosphere where the MW emission is 
formed can also contribute to this viewing angle difference. 
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